Exact trace formulae for a class of 

one-dimensional ray-splitting systems 

o" 



3 



Y. Dabaghian, R. V. Jensen and R. Bliimel 

Department of Physics, Wesleyan University, 

Middletown, CT 06459-0155 

February 1, 2008 

o\ 
o 
t> 

,__l , Based on quantum graph theory we establish that the ray-splitting 

trace formula proposed by Couchman et al. (Phys. Rev. A 46, 6193 

(^ ■ (1992)) is exact for a class of one-dimensional ray-splitting systems. 

^h, Important applications in combinatorics are suggested, 

-i— > 

§ ■ PACS: 05.45. +b,03.65.Sq,72.15.Rn 



>< 



Abstract 



1 Introduction 



Gutzwiller's trace formula, established in the late 1960s and early 1970s states 
that it is possible to obtain the level density of a bounded Hamiltonian dy- 
namical system with semiclassical accuracy, based entirely on the information 
provided by its classical periodic orbits [[[], ||]. According to Gutzwiller the 
density of energy states can be written as a sum over prime periodic orbits 
and their repetitions, 



1 oo 

p (£) = £ 6{E - E n ) = p(E) + —Re £ T p (E) £ A pv {E)e 

n ^ p v=\ 



^+*p(*0 



(1) 



where p(E) is the average density of states, S p is the classical action of the 
prime periodic orbit p, 

_ dS p (E) 

lp ~ dE [Z) 

is its period and (p p (E) is its Maslov phase. Gutzwiller derived the pre- 
exponential factors A pv in semiclassical approximation, expressing them in 
terms of the stability properties of the corresponding periodic orbits. 

The orbits used to construct the sum ([[[) are obtained at a given value of 
the energy E. On the other hand it is known that as the energy of a generic, 
nonhyperbolic system changes, the structure of the phase space changes and 
with it the set of periodic orbits. This phenomenon is called "phase-space 
metamorphosis" 0. Phase-space metamorphosis, in general, is accompanied 
by the creation and destruction of periodic orbits giving rise to the interesting 
phenomenon of ghost orbits Q. Therefore, in general, the sum in (|]) will 
change as a function of E in the sense that it may acquire or lose certain 
terms. However, apart from hyperbolic systems||], there exists an interesting 
class of systems, which are free of such metamorphoses of the phase space. 
These are the scaling systems, for which the action functional S P (E) for 
any periodic orbit p decouples into the product of an energy-dependent part 
f(E) and the "reduced action" S®, which depends only on the geometry of 
the periodic orbit, 

S P (E) = f(E)S° p . (3) 

Examples of such systems are various billiard systems, and also (with suitable 
definitions of scaling parameters) the hydrogen atom in strong electric and/or 
magnetic fields 0, |6|, |7[. For such systems neither the geometry of the phase 
space nor the geometry of the set of periodic orbits change with energy. 
Therefore the structure of the sum ([]]) can be defined once and is valid for 
all values of E. In such cases it is interesting to investigate the relationship 
between the fixed set of periodic orbits (the periodic orbit spectrum) and the 
quantum energy spectrum. 

Direct derivation of Gutzwiller's trace formula, as presented originally by 
Gutzwiller, is based on the saddle point approximation. This implies that this 
formula is meant to work only semiclassically, i.e. to predict only the highly 
excited energy levels with semiclassical accuracy. Indeed, for generic billiard 
domains, Gutzwiller's formula is not exact |1|. However, in certain special 
cases (|l]) is known to predict the entire energy spectrum exactly. An example 



is the harmonic oscillator. Moreover, exact "Gutzwiller-like" trace formulae 
do exist. A theorem by Anderson and Melrose || states that for any billiard 
there exists a set of pre-exponential factors A a , which makes the relationship 
([]]) exact. Other exact Gutzwiller-like trace formulae are obtained in the 
context of quantum graph theory |9|, [T(| ll| . We show that ray-splitting 



systems |T2[ provide further examples of exact Gutzwiller-like trace formulae. 
A specific example, a scaling one-dimensional step billiard, closely related 
to quantum graphs, is presented in the following section together with its 
generalized ray-splitting Gutzwiller formula. In Section III we present a 
proof for the exactness of the ray-splitting Gutzwiller formula. In Section 
IV we define the ray-splitting zeta function and relate it to cycle expansion 
techniques. In Section V we make use of the exactness of the ray-splitting 
Gutzwiller formula to prove a nontrivial combinatorial identity. In Section 
VI we discuss our results and conclude the paper. 

2 The model 

In this section we study the spectrum and the generalized ray-splitting Gutzwiller 



formula for the one-dimensional scaling step billiard ||12|, |13| , [14] , |15| jig] (Fig. 1) 



V(x) = { °' for < x < 6, ( 

W -\v = \E, for b < x < 1, y ' 

where A is the scaling constant and E is the energy of the system. In this 
paper we focus on the case A < 1. First results on this model were pre- 



sented in [jnj in the context of generalized Poisson formulae in ray-splitting 
systems. Despite its formal simplicity (H) can be used to illustrate many phys- 
ical and mathematical methods and ideas connected with the ray-splitting 



approach [|T2|, [13], |T^, [T5|, IBfl . We work in units such that ft — 1, the width of 
the potential well is 1 and the mass of the quantum particle is 1/2. Defining 
k = \/E and k = (3k, where 



p = VT^\, (5) 

it is elementary to obtain the exact quantum mechanical equation 

cos(kb) sink(l — 6)1 + — sin(A;6) cosk(l — 6)1 = (6) 

k 



for the energy levels E n of the system. They are determined by the roots k n 
of (H) according to E n = k\. 

It is more convenient to write (|6|) in the form 

sin(fcu;i) — rsm(kuj2) = (7) 

with 

U! = h + h, co 2 = h- h, h = b, l 2 = /3(1 - b) (8) 

and the reflection coefficient 

1-/3 



1 + P' 



(9) 



In general the two frequencies ui\ and uo 2 in (|7|) are not rationally related. 
Therefore (0) connects the physical problem of a scaling step potential with 
the mathematical theory of almost periodic functions ||17||. This means that 
all our exact results on the spectrum of the scaling step potential (|j) can be 
interpreted as theorems on the roots of doubly periodic functions. 

Solving (|7|) numerically, it is easy to obtain a large number of roots for 
studying statistical properties of the quantum energy levels as well as the 
relationships to the classical periodic orbits of the system. The latter goal is 
achieved by computing the Fourier image of the density of states defined as 

oo 

F(s) = Yl e~ lsk] - (10) 

3=1 

According to the Gutzwiller trace formula (|1|) , the Fourier transform provides 
a convenient tool for studying the orbit spectra of dynamical systems, since 
it produces pronounced peaks at those values of s that correspond to the 
actions of classical periodic orbits. In the case of the scaling system (f|), the 
actions S p in ([!]) are proportional to k, 



S P (E) = / k(x)dx = S"-k, (11) 

Jp 

and hence one expects the Fourier transform (10) of ([!]) to produce a <5-peak 



at s pv = vS)!, for every primitive periodic orbit p and its repetitions v. 



The result of the numerical evaluation of the sum ([10]) for this system 



is presented in Fig. 2. It shows a large number of narrow peaks. Most of 



them do not correspond to the standard (Newtonian) periodic orbits. This 
is immediately clear since in the case of the potential @ there exists only a 
single primitive Newtonian periodic orbit at any value of the energy above 
the potential step (see Fig. 1). 

The extra peaks in Fig. 2 are due to non-Newtonian periodic orbits fll2|, [13], 
[14], [15], [16|. They correspond to the non-Newtonian reflections off the sharp 
ray-splitting step. Together with the Newtonian orbits the non-Newtonian 
orbits account for every single peak in Fig. 2 for arbitrary values of the 
parameters A and b. Numerical computations indicate that the maxima of 
F(s) converge to <5-peaks in the limit when the number of roots included in 



the sum (|T0|) tends to infinity. This, in turn, suggests that there exists an 
exact formula of the type ([I]). That this is indeed the case is proved in Sect. 
Ill below. 

A natural generalization of Gutzwiller's trace formula, which includes 
the contributions from the non-Newtonian ray-splitting orbits, was obtained 
previously in fll2| . Speaking in terms of the step-potential (f|), instead of just 
a single orbit bouncing between x = and x — 1, a generic orbit may now 
be reflected off or transmitted through the ray-splitting boundary at x = b 
any number of times in arbitrary sequence. As a result, the set of primitive 
non- Newtonian orbits becomes infinite. In a one- dimensional system the 
numbers of these reflections and transmissions are the only characteristics 
of the orbits, and therefore any orbit can be characterized uniquely and 
completely by a binary sequence of symbols C and TZ that keep track of each 
reflection off the left (£) or the right (TV) wall of the potential well. The 
corresponding generalized Gutzwiller sum includes all primitive Newtonian 
and non-Newtonian periodic orbits and their repetitions, 



f.2T(p) r cr(p) 



V e iuS * (12) 



p(£) = p(£) + -Re£T p £[(-l)^)t 

71 v v=X 

where p = (b + /3(1 — b))/(2nk) is the average level density, r, defined by 
(H), is the quantum reflection coefficient, t = \/I — r 2 is the transmission 
coefficient, a(p) and 2r(p) are the number of reflections and transmissions of 
the primitive orbit p at the potential step, and xip) counts the total number 
of times the orbit reflects off the walls and off the potential step to the right 
of the ray-splitting boundary. Note that (— l)^( p ) = e llfp and explicitly defines 
the Maslov phase in ([!]). 



If we denote the actions of the shortest orbits (1Z and C) by S-ji and Sc 
respectively, the action S p of any orbit can be expressed as a sum 

S p = n c Sc + n n S n , (13) 

for certain integers ri£ and n-jt (generally different from a and r). The level 
density (|l|) contains only even powers of the transmission coefficient t, be- 
cause every periodic orbit transmits an even number of times through the 
ray-splitting boundary. 



3 Exactness of the ray-splitting trace formula 

Using quantum graph theory || [l(| |Tl| , it is possible to show that the ex- 
pression ([TJ) is exact. We prove this below after presenting some basic ideas 
of quantum graph theory. 

^From the perspective of quantum graph theory, the quantization of a 
particle in the potential (|4|) is treated as a scattering problem on the graph 

•l »2 »3 (14) 



with three vertices and two bonds described by the connectivity matrix || 

nana 

(° l ®\ 

C = 1 1 . (15) 

\o i o) 

On every bond connecting vertices i and j, one defines a free particle wave 
function ipij, which satisfies the following vertex conditions: 

4>iM = o) = &, 4>ij(£ = Lij) = <pj, (16) 

where £ is the coordinate along a particular bond of length L^, so that 
the wave functions on different bonds match on every vertex. The general 
solution satisfying the vertex conditions is 

^--fe-^-^^-***. (17, 



The presence of the coefficients fy allows us to generalize the formalism 
developed in pi, [TO, O. For the derivatives we have the continuity conditions 



£ C^'jM = U 3 ) = £ Ctf^tf = 0). (18) 

j<i j>i 

For the case of the potential (§), 

fa '« = sin[Wl-6)] Sin " ;/J(1 - 6 - fl '- (20) 

The matching and continuity conditions at vertex 2 result in 

/3tan(ife&)+tan[A;/3(l-&)]=0 (21) 

or 

sm[k(h + l 2 )] - r sin[k(h - l 2 )] = 0, (22) 

which is the same as ([!]). The lengths l\ and Z 2 , defined in (pj) turn out to be 
the weighted bond lengths. 

The same quantization condition can be obtained from considering the 
scattering process at every vertex of the graph (0). The vertex scattering 
matrix is given by 

a jiij' = djyCjiCif, (23) 

where the index i refers to the vertex under consideration, qfl are coefficients 
that depend on the physics of the scattering at the vertex i and CV,- are 
the matrix elements of the connectivity matrix determining the geometry of 
the graph. At the "dead end" vertices 1 and 3 of the graph (|14]) we have 
°jiii' = °j33j' = — l - ^ or ^ ne cen -tral vertex it is easy to show that 



The graph scattering matrix, describing the graph as a whole, is given by 

-1 
Da® 
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— (A) rT), <-) 



where the matrix 

M e«) < 26 > 

accounts for the phases accumulated along the bonds. The quantization 
condition m 1(| |H |, 

det (1 - S) = 0, (27) 

results in 



2i(hk+l 2 k) 



_ r jg«i* _ e 2ihk^ = L ( 2 g) 



This is the same as (|22|). 
With the help of 



OO 1 

lndet(l-S) = -]T ~Tr(S n ), (29) 



n=i n 



the quantization condition fl2~T|) can be written alternatively as a sum over the 
periodic orbits of the graph. Indeed, since the scattering matrix is defined 
geometrically using the graph connectivity matrix, its indices correspond to 
the vertices i and j connected by a bond if the matrix element Cy 7^ 0. The 
trace of the n-th power of this matrix is defined on the set of all the possible 
cyclical n-bond sequences. Using (p3|) and (|25| ) we obtain Tr(S 2n+1 ) = and 



Tr (S 2n ) 



\ ' /_ 1 )x r - t 2r e 2ifcL„ 5 ( 30 ) 



n c +n n =n 



where ric and n-ji give the number of times the left (£) and the right (TZ) 



bonds of flTID occur in the sequence, L n = n c h + n-nh-, & is the number 
of reflections from the middle vertex and It is the number of transmissions 
through it. Since the reflection coefficient coming from the scattering matrix 
(pip can be positive or negative, the factor (— l) x is needed to keep track of 
how many times it appears with a minus sign. It also keeps track of how 
many times a given orbit scatters off the walls. 

There are two possibilities for L n . Either it corresponds to a prime peri- 
odic orbit, or it corresponds to a bond sequence that retraces itself v times. 
In this case L n is v times the length of a single (primitive) traversal L p , 
and the pre-exponential factor is the z/-th power of the factor corresponding 
to the shortest closed bond sequence. Alternatively, these closed bond se- 
quences can be viewed as periodic orbits traced by a particle moving on the 
graph. 



On the other hand, from the analytical properties of det[l — S(E)] [O], 
it is easy to relate it to the spectral counting function N(E), 

1 1 °° 1 

N(E) = N W (E) - - + -Im £ -Tr (S n ) , (31) 

2 t n=1 n 

where N-w(E) is the average spectral staircase function (the Weyl term). 
Putting all these ingredients together, one arrives at the trace formula ([121) . 

The exactness of ( |i"2"D can be understood on the basis of the dramatic in- 
crease of the number of the primitive non-Newtonian orbits included in fll2|) . 
As the orbit length L p in ( |I2]) increases, non-Newtonian orbits proliferate 
exponentially, providing additional information about the structure of the 
potential V(x). 

It is interesting to compare ([12]) with the results produced by Gutzwiller's 
trace formula without ray splitting. For the step potential (^) there exits only 
one Newtonian periodic orbit with classical action Sn which bounces between 
the left and the right walls of the well (Fig. 1). Therefore fll2[) predicts 

r~n OO OO 

P (E)=p(E)+^J2 C0S (" S n)= T n E HS N -2nm), (32) 

v=\ m=—oo 

where 

S N = 2kb + 2kp{l -b) = S° N k, (33) 

T^ = 8Sn /dE and p = T^ /2n. Since only a single periodic orbit contributes, 
( |3"2"D predicts a periodic spectrum for a particle in the step potential (||). 
Moreover, the Fourier image 



r r i2 °° r n 

/ p(E) e~ lsk dE = [S N ] J2 5 [ mS °N ~ s \ ( 34 ) 



of the density of states produces <5-peaks at integer multiples of the reduced 
action S N of the (only) Newtonian orbit of the system. Figure 3 shows that 
the exact spectrum of the problem is not periodic, which illustrates that the 
trace formula (|I|) ((0), respectively) without non-Newtonian orbits predicts 
a wrong energy level distribution. 



A straightforward generalization of the ideas and procedures discussed 
above provides a proof of the exactness of (|l"2D for the whole class of iV-step 
scaled potentials 

V(x) = Vi = \E, b,^ <x<b u i = l,...,N, (35) 

where 60 — 0, b^ = 1 and A» is the scaling coefficient for the i-th interval 

[bi-i,bi\. 

4 Ray-splitting zeta function 



Despite its exponentially decreasing terms, ([12]) converges only conditionally 
due to the exponential proliferation of the non- Newtonian orbits. Therefore, 
one should specify a physically meaningful way of partial summation for this 
series. In practice, one could certainly consider the shortest periodic orbits 
in order to get an approximation for (|l|). Figure 3 presents the contribution 
of the 43 shortest periodic orbits. This corresponds to including all periodic 
orbits up to binary code length 7 (e.g. C1Z1ZCCCC). Figure 3 shows that 
the peaks give a very accurate representation of the actual positions of the 
roots. 



Since (|P2|) is a geometric series with respect to the repetition index u, this 



part of the summation can be performed immediately, yielding 

1 (_l)x(.P)+2T(p) r a(p) iS p (E) 

p(E)=p(E) + -ReY T v , w , „ „ , s -a ,^ - (36) 

Using the relation (g), the density of states can be written as 

p (E)=p(E)-hm-^lnZ(E), (37) 

where 

Z(E) = Y[ [l - (_i)x(p) t 2r( P ) rCT ( P ) e *s P (fc)j ( 38 ) 

p 

is an analog of the Fredholm determinant associated with the ray-splitting 

system ([|) and the sum (|T2"D , considered as a function of the coefficients 

r and t. One can consider a cycle expansion [18] of the product (38) in 



powers of r or t. A natural choice for the expansion variable would be the 
smaller one of r, t. Physically, this asymmetry determines whether reflection 
or transmission is the dominant process. 

10 



5 Combinatorics 

In quantum graph theory the representation of the quantum level density in 
the form of a Gutzwiller-like trace formula is exact || |H], |TIJ . In special cases 
both the level density and the trace formula can be evaluated analytically and 
give rise to combinatorial identities. This idea was successfully implemented 
by Schanz and Smilansky who obtained a host of new and nontrivial com- 
binatorial identities |l!|. Additional identities are generated whenever (|^) 



can be solved analytically. We illustrate the method by choosing f3 such that 
h = b and l 2 = (3(1 - b) = b in (H). In this case fl2|) becomes sin(2kb) = 0, 
solved by k n = mt/(2b). The corresponding periodic level density can be 
obtained directly using the conventional Poisson formula, 

It is interesting that the arguments of the exponents in (|39|) coincide with the 
actions of the repetitions of the Newtonian orbit, Sn = ^bk. This coincidence 
is due to the special choice of parameters, l\ = li = b assumed in this case. 



Alternatively, the level density (p9|) can be expressed via ([121) . Equating 



the pre-factors of terms with the same action results in the following sum 
rule 

1 = — V T p [(-l)X(p) t 2r( P ) r <x(p)j " _ ( 4 Q) 



pv 



Here T/v = 2b/k and the sum on the right-hand side is over all periodic 
orbits, Newtonian and non-Newtonian, that add up to the same multiple of 
the Newtonian action Spf. 

The sum rule ( |40|) can be recast into a combinatorial theorem on the set 
W\ of cyclically non-equivalent binary codes of even length A = 2M over the 
symbols £ and TZ in the following way. (i) For every word w G W\ compute 
the primitive time T w defined as T w = M/u w , where v w is the number of 
repetitions of the shortest sub-code in w. (ii) Scan each word w G W\ and 
assign 

w -> (-1)°- {r 2 f w (t 2 ) lw (41) 

according to the substitutions 

az ->t, nc^ t, cc ->r, nn^ -r. (42) 

11 



Then, with r + t = 1, we have 

£ T w (-l) a ™(r 2 Y w {t 2 y™ = 1, (43) 

w£W A 

which is equivalent to the sum rule (£|(]). Stated differently, (^) is the same 
as 

E (-IT-T W =( M ). (44) 

6 Discussion and conclusions 



Since Atle Selberg discovered his famous trace formula in 1947 [20], more 
exact trace formulae were found. There are many cases in which the geomet- 
rical information contained in the set of closed geodesies (periodic orbits) can 
be used to reconstruct the spectrum exactly Gutzwiller provided a physical 
theory which parallels these results. His formula points the way to establish 
approximate relationships, which involve physical rather than geometrical 
concepts. In particular, Gutzwiller uses the semi classical saddle-point ap- 
proximation, valid under certain physical conditions, in order to derive the 
pre-exponential factors in ([]]) in a form, which is valid for a wide class of 
dynamical systems. 

However, Gutzwiller's theory does not imply that these sums are neces- 
sarily approximate. There exist different approaches to establish exact rela- 
tionships between the spectra of operators and the spectra of periodic orbits. 
The above method, based on the intuitively clear idea of ray splitting, pro- 
vides an example which is physical, lies outside of the scope of Gutzwiller's 
approach, and is exact. 

This result immediately suggests many important applications. First, 
from the mathematical point of view, one can derive statements about the 
behavior of the zeros of a wide class of almost periodic functions. Second, 
the exactness of ([l2|) provides a convenient way to prove many combinatorial 
identities, expressible in terms of periodic orbits in (|j). Lastly, it provides 
a non-trivial way to obtain Feynman's path integrals in a well-defined limit. 
It is a natural idea to approximate arbitrary one- dimensional potentials by 
a step-like profile such as (|35l), for which ( |l"2"D is exact. Taking the limit in 



which the size of the steps tends to zero, one can approximate the shape of 

12 



any smooth potential with any accuracy. As the number of steps increases, 
the sum over the non-Newtonian orbits leads to Feynman's path integral. 

Y.D. and R.B. gratefully acknowledge financial support by NSF grants 
PHY-9900730 and PHY-9984075; Y.D. and R.J by NSF grant PHY-9900746. 
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Figure Captions 

Fig. 1: Step potential with non-Newtonian orbits C, CR.1Z and Newtonian 
orbit C1Z. 

Fig. 2: Fourier transform ( |i~0"D of the density of states of the step potential 
(fD with b = 0.7 and A = 1/2. About 10,000 states are included in the sum 
flTOl). Sharp peaks in the transform are located precisely at the actions of the 
Newtonian and non-Newtonian orbits. 

Fig. 3: Contribution of the 43 shortest Newtonian and non-Newtonian pe- 
riodic orbits (up to binary code length 7) to the density of states of the step 
potential shown in Fig. 1. The exact energy eigenvalues (+) are close to the 
locations of the peaks. For comparison the energy levels predicted by the 
Newtonian orbits alone are also shown (x). 
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